Challenges in generating density perturbations from a fluctuating inflaton coupling 
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We discuss the possibility of generating adiabatic density perturbations from spatial fluctuations 
in the inflaton decay rate which are due to quantum fluctuations of light moduli fields coupling to 
the inflaton. We point out that non-renormalizable operators, which lift the flatness of the moduli 
potential, play an important role for the density perturbations. In particular, the non-renormalizable 
terms give rise to a considerable damping of the fluctuations and thereby pose an obstruction to the 
construction of possible models. 



Inflation is the main contender for an explanation 
of the observed adiabatic density perturbations with a 
nearly scale invariant spectrum [1]. However, recently 
alternative mechanisms for generating the density pertur- 
bations have also been much discussed. In the curvaton 
scenario, iscocurvature perturbations of some light "cur- 
vaton" field are converted into adiabatic perturbations 
in the post-inflationary universe [2-4] . Another interest- 
ing proposal is that the perturbations could be generated 
from the fluctuations of the inflaton coupling to the Stan- 
dard Model degrees of freedom [5,6]. It has been argued 
that the inflaton coupling strength to ordinary matter, 
instead of being a constant, could depend on the vac- 
uum expectation values (VEV) of the various fields in 
the theory. If these fields are light during inflation their 
quantum fluctuations will lead to spatial fluctuations in 
the inflaton coupling strength. As a consequence, when 
the inflaton decays, adiabatic density perturbations will 
be created because fluctuations in the inflaton coupling 
translate into fluctuations in the reheating temperature. 

A particularly interesting implementation of this sce- 
nario is to consider the Minimal Supersymmetric Stan- 
dard Model (MSSM) plus an inflaton field [5]. There 
are many flat directions in the MSSM [7]. The moduli 
fields parameterizing these flat directions are light, and 
their quantum fluctuations during inflation produce fluc- 
tuations in the inflaton coupling. The degeneracy of the 
MSSM scalar potential is lifted by supersymmetry break- 
ing effects and by non-renormalizable operators. In this 
note we analyze the effect of these non-renormalizable 
operators on the produced density perturbations. We 
restrict the discussion to MSSM flat directions, but the 
results can easily be adopted to more general models. 

Assuming the inflaton is a gauge singlet, it can decay 
to normal matter through both normalizable and non- 
renormalizable interactions in the superpotential [5] 
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where <f> is the inflaton, h and h are the two Higgs dou- 
blets, and q and q c are quark and lepton superfields 
and their anti-particles. M is some cut-off scale which 



could be the GUT scale or the Planck scale. The ef- 
fective coupling for inflaton decay into Higgs fields is 
X h = A (l + jj + ...), where S is the VEV of one of 
the flat direction fields in the theory. The effective cou- 
pling for the inflaton decay to quarks is X q = jj with 
S = (q), (q c ), or (h). Effective couplings of this form can 
result from integrating out heavy particles. 

If a single decay channel dominates, the density con- 
trast is 6p/p ~ SX/X. 1 For X q ^> Xh the non- 
renormalizable couplings dominate, and the inflaton de- 
cays predominantly into quarks and anti-quarks. We will 
refer to this decay as indirect decay. This yields a density 
contrast 
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For light fields such that m| <C H 2 the quantum fluctu- 
ations are set by the Hubble scale SS ~ H. The density 
contrast required to explain the observed temperature 
anisotropy in the cosmic microwave background (CMB) 
radiation is Sp/p ~ 10~ 5 , which can be obtained for 
S 1 * ~ 1Q 5 H* provided there is no later damping of the 
fluctuations. Here and in the following, the subscript * 
denotes the corresponding quantity evaluated at the time 
observable scales leave the horizon, which is some 60 e- 
folds before the end of inflation. Since S* ^> H* the 
perturbation spectrum will be Gaussian. 

In the opposite limit A^ 3> X q the inflaton decays 
mainly into Higgs fields. We will call this the direct 
decay channel, as decay is mediated by renormalizable 
operators. The corresponding density perturbation is 
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x It is also possible that one channel is responsible for the 
fluctuations, but another is the main decay mode. Since the 
inflaton decay is exponential oc exp(— V / H) it would require 
some fine-tuning to obtain sizable density perturbations. 
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For this decay channel the density perturbations are in- 
dependent of S. The fluctuations will be Gaussian for 
S* 3> -H* and non- Gaussian in the opposite limit. It is 
important to remember that in Eqs. (2,3) the right hand 
side is always evaluated at the time of inflaton decay. 

As was pointed out in [5] the flat direction S must 
have a mass smaller than H* in order to obtain adi- 
abatic, nearly scale invariant, fluctuations. Therefore 
the Hubble-induced supergravity correction to the mod- 
uli mass should not exceed mj ~ 10 _1 i7» in order to 
have a successful scenario. This can be realized e.g. 
in no-scale supergravity models [8] with an Heisenberg 
symmetry imposed on the chiral fields in Kahler func- 
tion [9], or in models where inflation is driven by a D- 
tcrm. Apart from the Hubble induced mass terms there 
are also soft contributions from low energy supersym- 
metry breaking. However such contributions are small 
~ O(TeV), and do not pose any threat to the scenarios 
discussed in [5,6]. Finally, the flat directions are lifted 
by non-renormalizable operators in the supcrpotcntial of 
the form W = KS n /nM£~ 3 [10,11]. Within the MSSM 
with i?-parity conservation most of the flat directions are 
lifted by n = 4,5,6 non-renormalizable operators [11]. 
The flattest one is lifted by n — 9. If dominant, these 
non-renormalizable contributions can play an important 
role during and after inflation. 

Before discussing the effects of the non-renormalizable 
operators, let us start by analyzing the parameter space 
where they are sub-dominant and therefore can be ne- 
glected. The scalar potential for a flat direction can be 
written as 



V(S) - \m 2 s S 2 
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where k ~ 0(1), and mj ~ O(TeV) is the soft mass for 
the flat direction. We have assumed here that there is no 
Hubble induced mass correction during inflation. Such 
a mass term can be included, but it will not change our 
conclusions in any essential way. Requiring that the mass 
term dominates puts an upper bound on the VEV during 
inflation 



S* < (mM n - 3 /«) 1/( "~ 2) - 
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For the indirect decay channel S* ~ 10 5 iT», as fol- 
lows from Eq. (2). This translates into an upper 
bound on the Hubble constant during inflation: < 
10 6 , 10 8 , 10 9 GeV for respectively n = 4, 5, 6; here we 
have assumed M ~ M p \, k ~ 1 and ms ~ TeV. How- 
ever, the bound can be made stronger. Quantum fluc- 
tuations during inflation grow until a saturation value 
(S 2 ) w 3H 4 /8ir 2 m 2 s [12]. If one assumes that the VEV 
of S has a typical value 5* ~ \J (S 2 ), the bound on the 
Hubble constant becomes ff* < 10 3 , 10 8 , 10 10 GeV for 
n = 4, 5, 6 respectively. This result is independent of 
ms, except that ms ~ lO" 5 -/?* in order to get density 
perturbations of the observed size Sp/p ~ 10~ 5 . Tak- 
ing ms ~ TeV, this result is inconsistent for n = 4, i.e., 



domination of the mass term and density fluctuations of 
the observed size is incompatible. The results are only 
marginally consistent for n = 5. 

The density contrast generated through the direct de- 
cay channel is independent of see Eq. (3). A Gaus- 
sian perturbation spectrum requires < S*, with S*» 
bounded by Eq. (5) if the mass term is to dominate. 

We will now consider the opposite limit, in which the 
VEV of S is large and the non-renormalizable terms in 
the potential dominate. We should point out that the dy- 
namics of the potential Eq. (4) has already been studied 
in [2,3] in the context of MSSM curvaton models. A sim- 
ple analysis shows that during inflation the flat direction 
field condensate is slow-rolling in the non-renormalizable 
potential V NR ~ k 2 S 2( - n ~^ / Mp {n ~ 3) . The Hubble pa- 
rameter and the amplitude of the field can respectively 
be estimated as [3] 



S 1 * ~ k~^6^ M , 
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where S = 5S/S* - 

Now let us turn our attention to the density perturba- 
tions. The equations of motion for the homogeneous and 
the fluctuation parts are given by 



S + 3HS+ V'(S) =0, 
5S k + 3H6S k +(% + V"(S) )5S k = Q 
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where prime denotes derivative w.r.t. S. Since we are 
only interested in the long wavelength mode (k — > 0), 
using the slow roll approximation during inflation we get 



3HS + V'(S) =0, 
3H6S + V"(S)SS = 0. 



(10) 
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where we have omitted the subscript k, understanding 
that 5S is for the superhorizon mode. The evolution of 
the ratio of the fluctuations to the homogeneous mode of 
S in the non-renormalizable potential Vnr is 
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where i denotes the initial value. During inflation the 
zero mode obeys Eq. (10), which can be integrated to 
yield 



1 V"(Q \ \ -!/2(«-2) 

1 + 3P^3>T^H ' (13 » 



with AN the number of e-foldings after observable scales 
leave the horizon, and -ff e nd the Hubble constant at the 
end of inflation. In a slow roll regime V"(S*) < H 2 . For 
AJV = 60 we find S/Si - 0.7 for n = 4. We conclude 
that during inflation the damping of the perturbations is 
ncgligblc small. 
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After inflation the moduli field S slow-rolls (although 
marginally) with V" <~ H 2 (t) 1 and we can still use the 
slow roll approximation Eqs. (10, ff). The perturbations 
are further damped according to Eq. (12) until H ~ 
and the field starts oscillating in the quadratic potential. 
There is no further damping during this epoch of oscil- 
lations. The adiabatic density perturbations are gener- 
ated when the inflaton field decays, which happens when 
H ~ T0 ~ A 2 m0, with the mass of the inflaton. 2 
Thus, damping occurs between the end of inflation and, 
depending on which event happens first, inflaton decay 
or the onset of moduli oscillations. The total damping 
factor is 

D _ (^5/5')dccay ^ / max[r^,TOs] \ 2 

(5S/S) en d V -ffend / 

- max [A^ ocay , {ms/m^f] , (14) 

where we have used that at the end of slow roll in- 
flation -ff on d — m^. For the indirect decay channel, 
(A q )dccay = Sdccay/M. The correct level of density per- 
turbations is obtained if DH*/S* <~ 10 -5 , as follows 
from Eq. (2). This implies the limit Sdccay ~ 0.1 Af or 
> 10~ 2 to^. For the direct decay channel on the other 
hand (A^)dccay = A and density perturbations are of 
the observed magnitude provided DH*/M ~ 10 -5 , see 
Eq. (3). Assuming S* > H * , which assures Gaussian 
fluctuations, this requires Ao > 0.1 or ms ^ lO~ 2 m0. 

To conclude, the inclusion of non-renormalizable oper- 
ators in the potential of the moduli leads rather gener- 
ically to a considerable damping of the perturbations. 
This puts severe constraints on the parameters of the 
model. In particular, either the inflaton coupling to nor- 
mal matter should be rather large, A > 0.1, or the moduli 
mass should be large, ms > lO~ 2 m0. Another possibility 
is to consider low scale inflation. In this case the VEV 
of the moduli is small, and the non-renormalizable op- 
erators in the potential are subdominant. However, for 
the indirect decay channel we find that this appears to 
be inconsistent if the potential is lifted by n — 4 oper- 
ators, and only marginally consistent for n = 5. The 
constraints are milder for the direct decay channel. 
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2 We assume that the inflaton decays perturbatively. The 
results change if preheating occurs, and the inflaton decays 
through non-perturbative processes. 
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